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TITLE OF THE INVENTION 

METHOD AND APPARATUS FOR DETECTING 
PHYCOCYANIN-PIGMENTED ALGAE AND BACTERIA 
FROM REFLECTED LIGHT 

Inventors: Robert Vincent 

STATEMENT REGARDING GOVERNMENTAL INTEREST 

[0001] The present invention was made through funding from grant 

number NAG3-2629 from the National Aeronautical and Space Administration 
(NASA) to through the Ohio Aerospace Institute (OAI) as fiduciary agent. The 
United States Government has certain rights in the invention. 

BACKGROUND OF THE INVENTION 
[0002] The present invention relates to a method of detecting 

phycocyanin-pigmented algae or bacteria (cyanobacteria) in water from reflected 
light. These organisms are generally referred to as phycocyanins and 
cyanobacteria. 

[0003] In many instances it is desirable to be able to detect the 

presence microorganisms in water, particularly bodies of water that serve as a 
source for drinking water or that may serve as a site for recreation, such as for 
swimming boating, water sports and fishing. Many of these organisms in high 
concentrations can be harmful to the public and to the environment generally. 
[0004] It is particularly desirable to be able to be able to detect the 

presence microorganisms in water in a manner that is convenient and provides 
relatively immediate results so that the public may be warned or other actions 
taken to avoid or eliminate contamination of the assayed water. 
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[0005] The Laurentian Great Lakes have experienced toxin-producing 

blooms of the cyanobacterium Microcystis sp. on a number of occasions over the 
past decade, including a massive bloom in Lake Erie in 1995 that caused a 
variety of water quality problems and attracted broad public concerns (Taylor, 
1997; Brittain et al, 2000; Budd et al, 2002). 

[0006] Most freshwater systems in the world are affected by 

anthropogenic eutrophication, leading to undesirable increases in planktonic and 
benthic biomass. These phenomena often show large local differences and 
interactions with patterns of water flow. Among various problems, the amount 
and distribution of nuisance-forming cyanobacteria is of primary concern for 
water management. Cyanobacterial blooms may cause a variety of water quality 
problems, including dissolved oxygen depletion and subsequent fish kills, 
aesthetic nuisances (e.g., odors, scums, fish tainting, unsightliness), and 
unpalatable and possibly unsafe drinking water (Carmichael, 2001). Such 
problems can severely limit aquatic habitat, recreational activities, fisheries, and 
use of a water body as a potable water resource. Microcystis spp., a common 
bloom-forming species of cyanobacteria, was regularly documented in Lake Erie 
several decades ago, when the lake was heavily eutrophied as a result of 
anthropogenic activities (Makarewicz, 1993). Subsequent phosphorus 
abatement strategies initiated as part of the Great Lakes Water Quality 
Agreement have been largely successful resulting in a reduction in algal biomass 
and greater lake transparency. Despite these actions, blooms of Microcystis spp. 
have returned to Lake Erie, recurring each summer since 1995. The return of the 



blooms appears to coincide with the spread of invasive zebra mussels 
throughout Lake Erie and is possibly related to selective filtration of other 
phytoplankton by the mussels and rejection of Microcystis spp (Vanderploeg et 
al. 2001). During September, 1995, Lake Erie experienced a Microcystis spp. 
bloom resembling a thick slick of grass-green paint that extended over the entire 
surface of the western basin (Taylor, 1997; Brittain et al, 2000; Budd et al, 2002). 
Another notable bloom was reported in September, 1998 (Lake Erie LaMP, 
2000). These blooms are of special concern because at least some Lake Erie 
strains of Microcystis spp. produce the peptide hepatotoxin microcystin, which is 
harmful to waterfowl or other animals that might drink the untreated water 
(Brittain et al, 2000). Microcystin has also been identified as a tumor promoter, 
making long-term ingestion of even low levels of the toxin of concern (Falconer 
and Humpage, 1996; Carmichael, 2001). 

[0007] It would be of economic and public health value to be able to 

detect early stage (emergent) blooms of cyanobacteria, and Microcystis spp. in 
particular, especially if it is on a sufficiently timely basis for municipalities and 
recreation facilities to implement a response plan. It has been shown that remote 
sensing technology can be used to estimate the concentration and distribution of 
cyanobacteria through measuring the concentration of the pigment phycocyanin 
(Dekker, 1993), which is indicative of the presence of cyanobacteria. In waters 
off the southeastern coastal U.S. and the Gulf of Mexico, Subramaniam et al 
(2001) have applied a multispectral classification algorithm that employs data 
from the Sea-viewing Wide Field-of-View Sensor (SeaWiFS) for mapping blooms 
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of Trichodesmium spp., a marine cyanobacterium. In cyanobacteria, 
phycobiliproteins constitute the major photosynthetic accessory pigments 
(MacColl and Guard-Friar, 1987). Whereas in marine species the pink-colored 
phycoerythrin is the dominant accessory pigment, in fresh water taxa, such as 
Microcystis spp., phycocyanin is the dominant pigment (MacColl and Guard- 
Friar, 1987). With the availability of LANDSAT 7 imagery for every 16-day 
overpass period, the present invention is intended to develop a LANDSAT TM 
algorithm for detecting various levels of phycocyanin in western Lake Erie. The 
present invention also allows for the mapping of turbidity in Lake Erie and its 
tributaries, to investigate relationships between phycocyanin and turbidity. 
[0008] With the availability of LANDSAT Thematic Mapper (TM) 

imagery featuring an overpass cycle of every 16 days (8-day intervals, if both 
LANDSAT 5 and 7 are employed), one goal of the present invention was to 
develop a set of algorithms, methods of their use and devices for detecting 
cyanobacterial blooms in Lake Erie, based on a unique spectral signature 
produced by phycocyanin, a light-harvesting pigment complex ubiquitous among 
cyanophytes. 

[0009] In addition to the features mentioned above, objects and 

advantages of the present invention will be readily apparent upon a reading of 
the following description and through practice of the present invention. 

SUMMARY OF THE INVENTION 
[0010] In general terms, the present invention includes a method of 

determining the presence of phycocyanin algae or bacteria in water as well as a 
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measurement method followed by transmission of data to a remote processing 
site. 

[0011] The invention includes a method of determining the presence 

of phycocyanin algae or bacteria in water from light reflected therefrom. The 
method comprises the steps of: (a) obtaining a measurement of reflected light 
from the water, the measurement comprising a measurement of the respective 
amount of light in at least two, preferably five frequency ranges; and (b) relating 
the approximate amount of the phycocyanin in the water to the respective 
amounts of light by applying an algorithm relating the respective amounts of light 
in the at least two, preferably five frequency ranges to the amount of phycocyanin 
algae or bacteria in the water. This may be expressed in colonies per milliliter or 
otherwise through appropriate adjustment of the magnitude and dimensions of 
the algorithms described herein or generated by the present method. It will be 
understood that the expression of the amount of phycocyanins in terms of 
colonies per ml water is only one of several ways to express the amount, and 
that reference to mathematical equivalents refers to any mathematically or 
logically related algorithms or expressions. 

[0012] In a preferred embodiment, the measurement of reflected light 

from the water comprises the measurement of the respective amount of light in at 
least five frequency ranges: (i) from about 0.45 pm to about 0.52 pm (ii) from 
about 0.63 pm to about 0.69 pm; (iii) from about 0.76 pm to about 0.90 pm; (iv) 
from about 1 .55 pm to about 1 .75 pm and (v) from about 2.08 pm to about 2.35 
pm; and (b) relating the approximate amount of the phycocyanin in the water to 
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the respective amounts of light by applying an algorithm relating the respective 
amounts of light in the at least five frequency ranges to the amount of 
phycocyanin algae or bacteria in the water. 

[001 3] Preferably, the measurement of the amount of light in the at 

least five frequency ranges comprises the measurement, respectively, of: (i) 
LANDSAT Thematic Mapper ("TM") band 1, (ii) LANDS AT TM band 3, (iii) 
LANDSAT TM band 4, (iv) LANDSAT TM band 5 and (v) LANDSAT TM band 7. 
The algorithm may be any algorithm selected from the group consisting of: X * Ki 
- K 2 x (R31) + K3 x (R41) - K4 x (R43) - K 5 x (R53) + Ke x (R73) - K 7 x (R74) 
and equivalents wherein: 



[0014] 


X is the approximate amount of phycocyanin algae or bacteria 


expressed in micrograms per liter; 


[0015] 


K1 is a value in the range of from about 30 to about 60; 


[0016] 


K2 is a value in the range of from about 5 to about 15; 


[0017] 


K3 is a value in the range of from about 20 to about 35; 


[0018] 


K4 is a value in the range of from about 100 to about 1 30; 


[0019] 


K5 is a value in the range of from about 3 to about 10; 


[0020] 


K6 is a value in the range of from about 30 to about 50; 


[0021] 


K 7 is a value in the range of from about 5 to about 20; 


[0022] 


R31 is the value of LANDSAT TM band 3 divided by LANDSAT 



TM band 1, after subtraction for atmospheric haze separately in each band; 
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4 * 



[0023] R41 is the value of LANDSAT TM band 4 divided by LANDS AT 

TM band 1, after subtraction for atmospheric haze separately in each band; 
[0024] R43 is the value of LANDSAT TM band 4 divided by LANDSAT 

TM band 3, after subtraction for atmospheric haze separately in each band; 
[0025] R53 is the value of LANDSAT TM band 5 divided by LANDSAT 

TM band 3, after subtraction for atmospheric haze separately in each band; 
[0026] R73 is the value of LANDSAT TM band 7 divided by LANDSAT 

TM band 3, after subtraction for atmospheric haze separately in each band; and 
[0027] R74 is the value of LANDSAT TM band 7 divided by LANDSAT 

TM band 4, after subtraction for atmospheric haze separately in each band. 
[0028] Preferably, the K values are as follows: 

[0029] Ki is a value in the range of from about 45 to about 50; 

[0030] K 2 is a value in the range of from about 7 to about 1 1 ; 

[0031] K3 is a value in the range of from about 25 to about 35; 

[0032] K4 is a value in the range of from about 1 1 0 to about 1 20; 

[0033] K5 is a value in the range of from about 5 to about 8; 

[0034] K6 is a value in the range of from about 35 to about 45; and 

[0035] K 7 is a value in the range of from about 1 0 to about 1 5. 

[0036] Most preferably the K values are as follows: 

[0037] K1 is a value in the range of from about 46 to about 48; 

[0038] K 2 is a value in the range of from about 8 to about 10; 

[0039] K3 is a value in the range of from about 27 to about 30; 

[0040] K4 is a value in the range of from about 1 1 5 to about 1 20; 
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[0041] Ks is a value in the range of from about 6 to about 8; 

[0042] K6 is a value in the range of from about 38 to about 43; and 

[0043] K 7 is a value in the range of from about 1 3 to about 1 5. 

[0044] The method according to the present invention is such that the 



calculated value of phycocyanin-pigmented algae or bacteria (X) correlates to the 
actual measured amount of the phycocyanin in the water by an adjusted square 
correlation value (i.e., R 2 adjusted) in excess of 60% and as high as in excess of 
70%. 

[0045] The present invention may additionally comprise the step of 

generating a report of the approximate amount of the phycocyanin-pigmented 
species in the water. This may be done using electronics adapted to digitize and 
process the data using an appropriate algorithm, such as that described herein. 
For instance, the report may include an estimate of the number of the 
phycocyanin colonies per ml in the water. 

[0046] The method of the present invention may also include the step 

of transmitting data relating to the approximate amount of the phycocyanin in the 
water to a site remote from the site where the measurement takes place. This 
may be done using any transmission method including land line or wireless 
transmission. This is also used advantageously where the reflected light is 
sensed remotely by aircraft, satellite, boat or buoy. Processing of the data may 
take place at the site of light uptake or may be carried out at a remote location 
after transmission of the raw data. The estimated phycocyanin report may be 
sent to public authorities, such as police departments, fire and rescue 
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departments or life guard services to warn swimmers, boaters, sportsman or the 
public at large that a given body of water, or portion thereof, likely contains 
elevated/dangerous levels of phycocyanin algae or bacteria. 
[0047] The invention also includes an apparatus for determining the 

presence of phycocyanin algae or bacteria in water from light reflected therefrom, 
the device comprising: (a) a measurement device adapted to measure reflected 
light from the water, the measurement comprising a measurement of the 
respective amount of light in at least five frequency ranges: (i) from about 0.45 
pm to about 0.52 pm; (ii) from about 0.63 pm to about 0.69 pm; (iii) from about 
0.76 pm to about 0.90 pm; (iv) from about 1 .55 pm to about 1 .75 pm; and (v) 
from about 2.08 pm to about 2.35 pm; and (b) a processor capable of relating 
the approximate amount of the phycocyanin in the water to the respective 
amounts of light by applying an algorithm relating the respective amounts of light 
in the at least five frequency ranges to the amount of phycocyanin algae or 
bacteria in the water. 

[0048] Most preferably, the at least five frequency ranges comprise, 

respectively: (i) LANDSAT TM band 1, (ii) LANDSAT TM band 3 and (iii) 
LANDSAT TM band 4, such as (iv) LANDSAT TM band 5, and (v) LANDSAT TM 
band 7 and wherein the algorithm is any algorithm selected from the group 
consisting of: X - Ki - K 2 x (R31 ) + K 3 x (R41 ) - K4 x (R43) - K5 x (R53) + Ke x 
(R73) - K 7 x (R74) and equivalents wherein: 
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[0049] 


X is the approximate amount of phycocyanin algae or bacteria 


expressed in micrograms per liter; 


[0050] 


Ki is a value in the range of from about 30 to about 60; 


[0051] 


K2 is a value in the range of from about 5 to about 15; 


[0052] 


K 3 is a value in the range of from about 20 to about 35; 


[0053] 


K4 is a value in the range of from about 100 to about 130; 


[0054] 


K5 is a value in the range of from about 3 to about 10; 


[0055] 


K6 is a value in the range of from about 30 to about 50; 


[0056] 


K7 is a value in the range of from about 5 to about 20; 


[0057] 


R31 is the value of LANDSAT TM band 3 divided by LANDSAT 



TM band 1, after subtraction for atmospheric haze separately in each band; 
[0058] R41 is the value of LANDSAT TM band 4 divided by LANDSAT 

TM band 1, after subtraction for atmospheric haze separately in each band; 
[0059] R43 is the value of LANDSAT TM band 4 divided by LANDSAT 

TM band 3, after subtraction for atmospheric haze separately in each band; 
[0060] R53 is the value of LANDSAT TM band 5 divided by LANDSAT 

TM band 3, after subtraction for atmospheric haze separately in each band; 
[0061] R73 is the value of LANDSAT TM band 7 divided by LANDSAT 

TM band 3, after subtraction for atmospheric haze separately in each band; and 
[0062] R74 is the value of LANDSAT TM band 7 divided by LANDSAT 

TM band 4, after subtraction for atmospheric haze separately in each band. 
[0063] Preferably, the K values are as follows: 

[0064] Ki is a value in the range of from about 45 to about 50; 
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[0065] 


K 2 is a value in the range of from about 7 to about 1 1 ; 


[0066] 


K3 is a value in the range of from about 25 to about 35; 


[0067] 


K4 is a value in the range of from about 1 10 to about 120; 


[0068] 


K5 is a value in the range of from about 5 to about 8; 


[0069] 


Kq is a value in the range of from about 35 to about 45; and 


[0070] 


K 7 is a value in the range of from about 10 to about 1 5. 


[0071] 


Most preferably the K values are as follows: 


[0072] 


K1 is a value in the range of from about 46 to about 48; 


[0073] 


K 2 is a value in the range of from about 8 to about 10; 


[0074] 


K 3 is a value in the range of from about 27 to about 30; 


[0075] 


K4 is a value in the range of from about 1 1 5 to about 120; 


[0076] 


K5 is a value in the range of from about 6 to about 8; 


[0077] 


K6 is a value in the range of from about 38 to about 43; and 


[0078] 


K 7 is a value in the range of from about 13 to about 1 5. 


[0079] 


It is preferred that the apparatus is capable of performing such 



that the calculated value of phycocyanin correlates to the actual measured 
amount of the phycocyanin in the water by a correlation value in excess of 60, 
and most preferably by a correlation value in excess of 70. 
[0080] The apparatus may additionally include a report generator 

adapted to generate a report of said approximate amount of said phycocyanin in 
the water. Such a report generator may be any device that is adapted to place 
the data into a tangible medium, such as a printer, CD burner, flash memory, 
magnetic storage media, etc. 
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[0081] The apparatus may additionally include a transmitter adapted 

to transmit data relating to the approximate amount of the phycocyanin in the 
water from the processor to a site remote from the site where the measurement 
takes place. Such a transmitter may include those adapted to send data such as 
through land line or wireless transmission, including telephone, internet, cell 
phone, radio and the like. 

[0082] The measurement device may be any device adapted to sense 

and record and/or transmit the light frequencies described above. Examples 
include photosensors, cameras, digital cameras and video cameras, etc. 
[0083] The processor may be any data processing device having 

programming instructions for applying the algorithm, such as preferably a 
microprocessor. 

[0084] It is preferred that the algorithm comprises a linear relationship 

between the approximate amount of the phycocyanin in the water and sum of (a) 
the ratio of the first frequency to the second frequency and (b) the ratio of the 
second frequency to the third frequency. 

[0085] The measurement device may be placed in any position from 

which it can sense the required light frequencies, such as on a buoy, a boat, a 
light house or similar dedicated tower structure, an elevated lifeguard house. 
The measurement device may also be in the form of a handheld device, such as 
a camera connected to a processor for processing the recorded light frequencies, 
the device may also be in the form of a device similar to a personal digital 
assistant with light recording and processing functions. 
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[0086] Another variation of the invention is a system using 

transmission of light measurement data to processor at a different location, 
recognizing that the processing may be done at a different location than the light 
sensing/recording. 

[0087] In general terms, this variation is a system for determining the 

presence of phycocyanin algae or bacteria in water from light reflected therefrom, 
the device comprising (a) a measurement device adapted to measure reflected 
light from the water, the measurement comprising a measurement of the 
respective amount of light in at least five frequency ranges: (i) from about 0.45 
pm to about 0.52 pm; (ii) from about 0.63 pm to about 0.69 pm; (iii) from about 
0.76 pm to about 0.90 pm; (iv)from about 1.55 pm to about 1.75 pm; and (v) 
from about 2.08 pm to about 2.35 pm; and (b) a processor at the remote site and 
capable of relating the approximate amount of the phycocyanin in the water to 
the respective amounts of light by applying an algorithm relating the respective 
amounts of light in the at least five frequency ranges to the amount of 
phycocyanin algae or bacteria in the water. 

[0088] The invention also includes a method of developing an 

apparatus for determining the presence of phycocyanin algae or bacteria in water 
from light reflected therefrom, the device comprising (a) obtaining a 
measurement of reflected light from the water, the measurement comprising a 
measurement of the respective amount of light of at least two frequencies; (b) 
developing an algorithm relating the respective amounts of light in the at least 
two frequencies to the amount of phycocyanin algae or bacteria in the water 
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through linear regression analysis; (c) producing a processor capable of relating 
the approximate amount of the phycocyanin in the water to the respective 
amounts of light by applying an algorithm relating the respective amounts of light 
in the at least five frequency ranges to the amount of phycocyanin algae or 
bacteria in the water; and (d) providing a measurement device adapted to 
measure reflected light from the water and adapted to provide data relating to the 
measurement to the processor. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0089] Novel features and advantages of the present invention, in 

addition to those mentioned above, will become apparent to those skilled in the 
art from a reading of the following detailed description in conjunction with the 
accompanying drawings summarized as follows: 

[0090] Figure 1 is a LANDSAT 7 ETM+ image (bands 1 , 2, and 3 

displayed as blue, green, and red, respectively) of part of the July 1 , 2000, Path 
20 Row 31 frame, showing collection sites of 30 water samples as black dots in 
Maumee Bay, Ohio. North is to toward the top, as taken in accordance with one 
embodiment of the present invention; 

[0091] Figure 2 is a graph of actual and predicted values of relative 

phycocyanin concentration (PC), according to the best model derived from dark- 
object-subtracted single bands of LANDSAT TM imagery for Lake Erie water 
samples collected on the same date as the LANDSAT 7 overpass (July 1, 2000), 
in accordance with one embodiment of the present invention; 
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[0092] Figure 3 is a graph of actual and predicted value of relative 

phycocyanin concentration (PC), according to the best model derived from dark- 
object-subtracted spectral ratios of LANDSAT TM data for Lake Erie water 
samples collected on the same date as the LANDSAT 7 overpass (July 1 , 2000), 
in accordance with one embodiment of the present invention; 
[0093] Figure 4 is a graph of actual (from the water samples collected 

on September 27, 2000) and predicted values of relative phycocyanin 
concentration (PC) showing application of the best LANDSAT 7, July 1, 2000 
spectral ratio model for relative phycocyanin concentration (PC) to the LANDSAT 
5 dataset (September 27, 2000), in accordance with one embodiment of the 
present invention; 

[0094] Figure 5 is a graph of actual (from the water samples collected 

on September 27, 2000) and predicted values of relative phycocyanin 
concentration (PC) showing application of the best LANDSAT 7, July 1, 2000 
single band model for relative phycocyanin concentration (PC) to the LANDSAT 
5 dataset (September 27, 2000), in accordance with one embodiment of the 
present invention; 

[0095] Figures 6 is a graph of relative phycocyanin content (PC) 

displayed as red (8.06-9.25 pg/l) to blue (0-5.17 pg/l), from the July 1, 2000 best 
spectral ratio model, applied to the July 1, 2000, LANDSAT 7 frame, in 
accordance with one embodiment of the present invention; 
[0096] Figure 7 is a photograph illustrating phycocyanin content (PC) 

displayed as red (10.31-15.77 Mg/I) to blue (0-2.50 |jg/l). from the July 1, 2000 
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best spectral ratio model, applied to the September 27, 2000, LANDSAT 5 frame, 
in accordance with one embodiment of the present invention; 
[0097] Figure 8 is a graph of percent reflectance (%) versus 

wavelength (nanometers) of absorption features of lake water containing 
primarily chlorophyll a (top curve) and water from another lake containing both 
phycocyanin and chlorophyll a (after Green, 2003). Wavelength limits of 
LANDSAT TM bands 1 (450-520 nm), 2 (520-600 nm), and 3 (630-690 nm) are 
shown as vertical bars, but limits for TM bands 4 (760-900 nm), 5 (1,550-1,750 
nm), and 7 (2,080-2,350 nm) are not shown. Spectral locations of absorption 
bands for chlorophyll a, phycocyanin, and dissolved organic matter 
(DOM)/chlorophylls are identified, in accordance with one embodiment of the 
present invention; 

[0098] Figures 9 is a graph of actual turbidity plotted versus actual PC 

for the July 1 , 2000 water samples, in accordance with one embodiment of the 
present invention; 

[0099] Figure 10 is a graph of actual vs. predicted turbidity from the 

July 1, 2000 turbidity model, applied to the July 1, 2000 LANDSAT 7 frame for 
P20R31, in accordance with one embodiment of the present invention; 
[00100] Figure 11 is a photograph illustrating turbidity (left) and PC 

(right) images of the Maumee River Mouth subregion (SW corner of Lake Erie) of 
the July 1, 2000 LANDSAT 7 ETM+ frame), in accordance with one embodiment 
of the present invention; 
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[00101] Figure 12 is a photograph illustrating phycocyanin content of 

P19R31 of Western Lake Erie (red is highest PC, from 4.98-12.00 pg/L) for 
LANDSAT 7 ETM data of July 16, 2002, according to the July 1, 2000 model, in 
accordance with one embodiment of the present invention; 
[00102] Figure 13 is a photograph illustrating phycocyanin content of 

P19R31 of Western Lake Erie (red is highest PC, from 4.98-12.00 pg/L, and dark 
blue is zero) for LANDSAT 7 ETM data of August 1, 2002, according to the July 
1, 2000 model, in accordance with one embodiment of the present invention; and 
[00103] Figure 14 is a photograph illustrating Phycocyanin content of 

P20R31 of westernmost Lake Erie (red is highest PC, from 11.57-12.51 pg/L, and 
dark blue is lowest PC, 0.00-5.62 pg/L) for LANDSAT 7 ETM data of August 8, 
2002, according to the July 1, 2000 model, in accordance with one embodiment 
of the present invention. 

DETAIL DESCRIPTION OF PREFERRED EMBODIMENT(S) 
[00104] The preferred system herein described is not intended to be 

exhaustive or to limit the invention to the precise forms disclosed. They are 
chosen and described to explain the principles of the invention and the 
application of the method to practical uses so that others skilled in the art may 
practice the invention. 

[00105] The present invention includes the use of algorithms developed 

from LANDSAT 7 ETM+ data for the July 1 , 2000 overpass and LANDSAT 5 TM 
data for the September 27, 2000 overpass for Path 20 Row 31 (including Toledo, 
Ohio) to measure relative phycocyanin content (PC) and turbidity in the western 
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basin of Lake Erie. Water samples were collected from discrete hydrographic 
stations arranged in a 20 km x 4 km grid adjacent to the Ohio shoreline during a 
6-hour period spanning each of the two LANDSAT overpasses. The samples 
were analyzed for Chlorophyll a content and turbidity. In addition, the 
concentration of phycocyanin, a light-harvesting pigment associated with 
cyanobacteria, was estimated from the ratio of phycocyanin :chl a in vivo 
fluorescence (IVPF/IVCF). A dark-object-subtracted, spectral ratio model derived 
from the July 1 , 2000 data was found to be the most robust, when applied to the 
September 27, 2000 data. The same July 1 , 2000 model (or algorithm) for PC 
was then applied to LANDSAT 7 ETM+ frames for July 16 and August 1, 2002 of 
the Path 19 Row 31 frame (including Cleveland, Ohio) and to the August 8, 2002 
frame of Path 20 Row 31. Moderate, very low, and high PC values were 
detected in the western basin of Lake Erie on July 16, August, 1, and August 8, 
2002, respectively. On September 17, 2002, local media reported a large 
Microcystis sp. bloom in the western basin. The high PC values on August 8, 
2002 may have been an early detection of the large Microcystis sp. bloom that 
was reported 5 weeks later. The PC algorithm derived in this study will improve 
our understanding of the temporal and spatial dynamics of cyanobacterial bloom 
formation in Lake Erie and other systems. It may also serve to alert 
municipalities to the presence of potentially toxic bloom events by sending 
relevant data and reporting these results. 

[00106] The invention also includes a system using an algorithm for 

converting LANDSAT TM multispectral signals into images showing different 
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values of phycocyanin colonies per milliliter of water. This system and method 
were tested in Lake Erie and its wider tributaries detect phycocyanin algae in the 
waters of Lake Erie to analyze the changes in water populations as they may 
affect human and wildlife activities. By gathering water samples during the 
period of time the satellite passes over Lake Erie and applying test kits, the level 
of phycocyanin populations was determined. The preferred algorithm combines 
2 ratios of three of the six spectral bands within silicon detector range (one to 
determine chlorophyll and the other turbidity). 

[00107] The method of the present invention may be carried out using 

any sensing appropriate light sensing devices adapted to capture the algorithm- 
relevant frequencies as described herein, including satellite and surface sensors 
(such as a wireless multispectral phycocyanin detector, or WIMCOD) for 
detection of phycocyanin, as well as for phycocyanin pigment, which is found in 
several types of toxic algae. Phycocyanin is a natural pigment that is not 
chlorophyll pigment, but which is found in some blue-green algae. 
[00108] An algorithm that may be used in the present invention, which 

may be carried out by computer instructions for producing a particular type of 
image that can be used to map a particular substance from a remote sensing 
platform in space, in an aircraft, or on the ground, may be determined as follows. 
[00109] LANDSAT Thematic Mapper (TM) is a sensor that has 8 

spectral bands, 6 of which have a 30-meter spatial resolution and which detect 
visible and infrared radiation (sunlight) reflected off the Earth's surface. The 
following bands were employed, with the wavelength limits (in micrometers, or 
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Mm) of their spectral band-widths given below for the LANDSAT 7 version of TM, 
called ETM+, and the LANDSAT 4 and 5 versions, called TM: 





P 1 f 






ElMg^^'&Fa'i|tta'5"d^J.dt^ 




-it 


Bandwidth (u) Full Width - Half Maximum 


Senso 
r 


Band 1 
Plot 
Data 


Band 2 
Plot 
Data 


Band 3 
Plot 
Data 


Band 4 
Plot 
Data 


Band 5 
Plot 
Data 


Band 6 
Plot 
Data 


Band 7 
Plot 
Data 


Band 8 
Plot 
Data 


TM 


0.45- 
0.52 


0.52- 
0.60 


0.63- 
0.69 


0.76- | 
0.90 


1.55- 
1.75 


10.4- 
12.5 


2.08- 
2.35 


N/A 


ETM+ 


0.45- 
0.52 


0.53- 
0.61 


0.63- 
0.69 


0.78- 
0.90 


1.55- 
1.75 


10.4- 
12.5 


2.09- 
2.35 


.52 - .90 



[00110] For instance, band 2 of the LANDSAT 7 version of the TM 

sensor (called ETM+) has wavelength limits of 0.53-0.61 Mm, band 3 has limits of 
0.63-0.69 and band 4 has limits of 0.78-0.90 pm. When mapping 
phycocyanin pigment, phycocyanin algae in Lake Erie and its tributaries with 
LANDSAT 7 data, it had to be determined which or how many of bands 1-5 and 7 
(which have 30-m spatial resolution and relatively narrow spectral bands, as 
opposed to the 60-m spatial resolution of band 6 and the relatively wide band- 
width of the 1 5-m-resolution band 8) to use. A mathematical procedure (multiple 
regressions) was applied to seek the best combinations of those bands to 
correlate with the target phycocyanin concentration. 
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[00111] It was determined that the use of the single band radiances 

(even if they were reduced to spectral reflectances from theoretical atmospheric 
models) as inputs to this procedure, the resulting algorithm would not be very 
robust (i.e., repeatable under different solar illumination and atmospheric 
conditions). Therefore, spectral ratios (ratios of spectral bands, after empirical 
correction for atmospheric haze through a process referred to as "dark object 
subtraction," were input to the mathematical procedure for each pixel from which 
a water sample had been collected. These 15 non-reciprocal ratios (R21, R31, 

R32, R41, R75) became the dependent variables and phycocyanin became 

the independent variable, which was the result of lab analysis of the water 
samples. For the LANDSAT 7 overpass, 30 water samples were collected, which 
were measured for both phycocyanin and phycocyanin content. The best 
subsets of spectral ratios were determined, and then the ones with the highest R 2 
(Adjusted) values were tested to see if they passed the Durbin-Watson test. The 
model with the highest R 2 (Adjusted) that also passed the Durbin-Watson test 
was the model that was considered to be the best. 

[001 1 2] The present invention thus takes advantage of two technology- 

based approaches: the use of remote sensing as a tool to study regional-scale 
aquatic ecosystem dynamics, and the use of sensitive fluorescence methods to 
identify and quantify algal pigments. The present invention has potential 
economic and public health value to municipalities located along the lake 
because it aims toward measuring the abundance of cyanobacteria on a 
sufficiently timely basis to allow municipalities and recreation facilities that 
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depend on Lake Erie for drinking water to respond to the threat. The best 
resulting algorithm was applied to Lake Erie and its tributaries, as well as to 
some small inland lakes in Northern Ohio, though only Lake Erie results are 
reported here. It also has potential for assisting future assessment and 
monitoring of cyanobacteria blooms, water quality, and aquatic ecosystem health 
in other regions and, perhaps, on a global scale. These results were obtained in 
a freshwater lake, however, and the algorithm may require changes for use in 
marine environments. 

[001 13] Powerful remote sensing techniques have become available in 

the last two decades that facilitate practice of the present invention to determine 
large-scale biological processes in difficult environments. At least four satellites 
have been commonly used for chlorophyll and phytoplankton mapping: AVHRR, 
SEAWIFS, and MODIS, all of which have spatial resolutions that range from 250- 
1,000 meters in pixel size. Where one is interested in results beyond the Great 
Lakes, to other fresh water lakes and their tributaries, as well as small inland 
lakes, the 30-meter resolution of the six visible/reflective IR spectral bands of 
LANDSAT TM and ETM+ are preferably selected. However, LANDSAT TM data 
has traditionally had one exceptional disadvantage: though data are collected by 
the LANDSAT satellites with a 16-day frequency (8-day frequency for two 
LANDSAT satellites), the data were not readily available to civilian scientists in 
less than approximately 60 days following the data collection. The availability of 
LANDSAT TM data within 24-48 hours through the OhioView consortium (a 
remote sensing consortium of eleven of Ohio's public research universities) 



23 



permits non-government scientists to perform time-sensitive research with 
LANDSAT data for the first time since ERTS I (later called LANDSAT I) was 
orbited in 1972. 

[001 14] An experiment was conducted on July 1 , 2000, a day of 

LANDSAT 7 overpass. Thirty sites with a grid spacing of 2 km over three, 20-km- 
long lines were surveyed (with GPS coordinates) in Maumee Bay, which is 
located at the southwestern corner of Lake Erie, with surface water samples 
collected at each site. Figure 1 shows a natural color image of Maumee Bay (a 
small part of the July 1 , 2000 LANDSAT ETM+ frame), with the location points of 
all 30 water samples collected within 3 hours of the overpass, about 30 minutes 
before noon, EDT. This followed the Presidential Order that permitted the un- 
fuzzing of the GPS satellites in the U.S. in May, 2000, resulting in precise 
location data correct to within approximately 3 meters in both Easting and 
Northing, which is much less than the size of a LANDSAT TM pixel (28.5 m x 
28.5 m) for the six reflective bands (1-5 and 7). In situ pH value, temperature, 
and turbidity were also measured at each site. Sampling was also conducted on 
September 27, 2000, coinciding with a LANDSAT 5 overpass, with 20 sites, but 
without turbidity data. The data collection sites were only approximately the 
same as for the first LANDSAT overpass, but GPS measurements were made at 
each site. 

[001 1 5] Chlorophyll a retained on a 0.2 j^m PCTE filter (GEOsmonics, 

Minnetonka, Minnesota) was determined by fluorometric analysis following 
extraction in 90% acetone at -20 °C for 24 hours (Welschmeyer 1994) using a 
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TD 700 fluorometer (Turner Designs, Sunnyvale, CA). Phycocyanin in vivo 
fluorescence (IVPF) was measured on samples using narrow-band interference 
filters (Andover Corp., NH) with excitation at 630 ± 10 nm and emission at 660 ± 
10 nm (Watras and Baker, 1988). In vivo Chlorophyll a fluorescence (IVCF) was 
measured on the same fluorometer, but using a filter set with excitation and 
emission wavelengths of 430 ± 10 nm and 680 ±10 nm. In the present 
invention, one may use the product of the ratio of IVPF/IVCF and extractive 
chlorophyll a to provide an estimate of the relative concentration of phycocyanin 
(PC = Relative Phycocyanin Content). PC values are proportional to absolute 
phycocyanin content. Cultures of Microcystis aeruginosa were used as 
standards to calibrate the fluorometer response for IVPF. Used in this capacity 
were M. aeruginosa UTCC 124, a non-toxic congener from the University of 
Toronto Culture Collection and M. aeruginosa strain LE-3, an isolate from Lake 
Erie collected during the 1995 bloom (Brittain et al, 2000). 
[001 16] LANDSAT 7 and LANDSAT 5 TM data were processed using 

the ERMAPPER image processing software package and Minitab statistical 
software package, both of which are commercially available. Dark object 
subtraction (Vincent, 1997), which is briefly described below, was applied to each 
band to reduce the effects of atmospheric haze from time one to time two. 

[001 1 7] The spectral radiance detected by the ith spectral band sensor 

can be approximated by the following equation: 
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[001 18] L, = \{sa\ p x + p\ )dX * (sa^ Pi + fit X'W " >W ) = ' «, Pt + A 

(Eqn. 1) 

[001 19] where Z f = the spectral radiance detected by the ith spectral 

band sensor for a given pixel on the Earth's surface. 

[001 20] k Uppgt = upper wavelength limit of ith spectral band. 

[00121] A Uwer = lower wavelength limit of ith spectral band. 

[00122] s = the shadow/slope factor that is 0 for total shadow and 1 for 

no shadow or slope in that pixel. 

[001 23] a. -a (X Vpper - X Lower ) =wavelength-dependent, multiplicative 

factor that includes atmospheric transmission and sensor gain averaged over the 
ith band. 

[00124] />.= spectral reflectance of the Earth's surface averaged over 

the ith band. 

[00125] p i = wavelength-dependent, additive factor that includes 

atmospheric 

[00126] haze averaged over the ith band and sensor additive offset. 

[00127] The approximation in Equation 1 is that the spectral band-width 

is much narrower than the width of most spectral absorption features of the 
object being observed. This approximation gets increasingly better with narrower 
spectral band-widths. All of the chemical composition information to be obtained 
from the surface of the Earth in a given pixel is found in the p. term, and in none 
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of the other terms. Thus, it is desirable to remove the additive term, which can 
be done simply by histogramming each spectral band and determining the 
minimum digital number (DN) found in all (or nearly all) of the pixels in the image. 
One less than this minimum DN is taken to be the value of the p. term (the dark 

object for the ith spectral band) in Equation 1 , which is the same for the entire 
frame of image data, under clear-sky conditions and a reasonably flat terrain or 
water surface. When one less than the minimum DN number is taken to be the 
dark object for a given spectral band, there will be no division by zero when that 
spectral band is employed as the denominator of a ratio. There is a different 
dark object for each spectral band, and the magnitude of the dark object is 
greatest for TM band 1 , less for band 2, even less for band 3, and sometimes as 
low as zero for band 4. Dark objects for TM bands 5 and 7 are usually zero. 
[00128] Once the dark object is found for the ith band, it is subtracted 

from all other pixels in the scene, yielding the equation for dark-object-subtracted 
radiance II. : 

[001 29] L. = L i - p. = SUM Eqn. 2 

[001 30] Although a i is approximately constant over the whole frame, 

under the assumed conditions of clear sky and small elevation differences 
(reasonably flat) throughout the image, both s and p t vary from pixel to pixel. 

Because s holds no information about the chemical composition of the pixel, it is 
highly desirable that it should be eliminated, which can be accomplished with a 
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spectral ratio (Vincent, 1997). A spectral ratio, denoted asi? fJ , is the quotient of 
the dark-object-subtracted radiances in the ith and jth spectral bands, or 

_ 4 sa iPi , a i\rPi\ 

[001 31] Rtj =-h = — = (—)(—) Eqn. 3 

L sa p a p 

J J^J J rj 

[00132] Note that the shadow/slope factor is cancelled out, and what is 

left are two multiplicative terms, the first of which is an atmospheric/multiplicative 
gain factor that is the same throughout the image, under the assumed conditions. 
The second term is a ratio of reflectances in the ith and jth bands for the pixel 
being observed and carries all the information about chemical composition. The 
spectral ratio of Equation 3 has been demonstrated many times on land to be 
more robust to solar illumination, atmospheric, and sensor parameter changes 
than any parameter based on linear combinations of single spectral bands 
because both p. and s have been removed by the dark-object-corrected ratio 
process. 

[00133] t. for six single bands of TM and the fifteen non-reciprocal 

spectral ratios (R u ) that can be produced from those 6 single bands were 

extracted from the image data for each of the pixels that contained the water 
sample collection sites for the July 1 , 2000 overpass of LANDSAT 7 (30 
samples) and the September 27, 2000 overpass of LANDSAT 5 (22 samples). 
Both types of data (single band combinations and spectral ratio combinations) 
were used to construct multiple regression models describing the relationship 
between the LANDSAT TM data and each of the measured values of turbidity, 
chlorophyll a and the ratio IVPF/IVCF, as well as the relative phycocyanin 
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concentration (PC), which was estimated for each water sample from the product 
of chlorophyll a content and the IVPF/IVCF ratio. 

[00134] Moving beyond any pre-conceived notions concerning the best 

single bands or the best spectral ratios for mapping the parameters of interest 
(PC and turbidity), the present invention may be arrived at by inputting some or 
all of them and allowing the multiple regression method sort the best ones out. 
The best multiple regression models for PC and turbidity were sought separately 
for linear combinations of single band and linear combinations of spectral ratios. 
It was our expectation, from Equation 3 and past experience on land 
experiments, that the spectral ratio models would be more robust (applicable to 
data collected from later dates) than the single band models. These 
expectations were proved correct, as shown in the next section. 
[00135] The mathematical models for both LANDSAT 7 and LANDSAT 

5 images were generated by using a step-wise linear regression method, with a 
Durbin-Watson statistical test for autocorrelation, which has been detailed by 
Vincent (2000). The Durbin-Watson (DW) test (Durbin and Watson, 1951) 
provided assurance that autocorrelation was not important in the parameters 
used for each final model reported. Although the DW test is more commonly 
employed in judging autocorrelation in multiple regression models of time series 
(multiple time inputs), it is equally useful forjudging autocorrelation in regression 
models from multiple spectral band inputs. The test applies to regression 
models produced from any type of multiple inputs. The best model of each type 
of input (single bands and spectral ratios) was applied to the image the model 
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came from, as well as to another image (with a different overpass date in year 
2000), to test how well and how robustly each type of regression model 
performed. 

[00136] The best PC and turbidity models were applied from July 1 , 

2000 to LANDSAT 7 data collected on July 16, 2002 and on August 2, 2002 of 
the adjacent frame area to the East (P. 19, R. 31). It was later reported by local 
media that a Microcystis spp. bloom occurred in the western basin of Lake Erie 
late in the summer of 2002. 

Results and Discussions 

[001 37] For the LANDSAT 7 image collected on July 1 , 2000, strong 

correlations were found for both the best single band model and the best spectral 
ratio model. The best single band model, which had an R 2 (Adjusted) =73.8% 
and a standard error of S = 0.6402 micrograms/liter (about 16% of the total PC 
range for the July 1 , 2000 overpass), is given by the following equation for 
relative phycocyanin content: 

[00138] PC=0.78 - 0.0539(B1) + 0.176(B3) - 0.21 6(B5) + 0.117(B7) 

Eqn. 4 

[00139] where B1, B3, B5, and B7 stand for dark-object-subtracted 

digital numbers of LANDSAT TM band 1, band 3, band 5 and band 7, 
respectively. 

[00140] In the preferred embodiment, the best spectral ratio model, 

which had an R 2 (Adjusted) = 77.6% and S = 0.5921 micrograms/liter (about 
14.8% of the total PC range for the July 1 , 2000 overpass), is given by 
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[00141] PC (MQ/L) = 47.7 - 9.21 (R31 ) + 29.7(R41) - 1 18(R43) - 

6.81 (R53) + 41.9(R73) -14.7(R74) 
Eqn. 5 

[00142] wherein RIJ (formerly called R tJ ) stands for the dark-object- 

subtracted spectral ratio of the Ith band over the Jth band. Only the spectral 
ratio model definitively passed the DW test; the single band model was 
undetermined in the DW test. Predicted values by the single band model and the 
spectral ratio model of phycocyanin for the July 1, 2000 overpass of LANDSAT 7 
are plotted in Figure 2 and Figure 3, respectively, versus PC for each of the data 
collection sites in Lake Erie. The same two models were then applied to the 
LANDSAT 5 data obtained on September 27, 2000. The spectral ratio model did 
reasonably well, as shown in Figure 4, and the predictive error of PC from the L7 
overpass model (July 1 , 2000) was 2.002 micrograms/liter for each value of PC 
measured in the L5 overpass (September 27, 2000) water sample data set, or 
about 18.2% of the total range in PC for the September 27, 2000 water samples. 
As shown in Figure 5, the single band model produced worse results than did the 
spectral ratio model, proving that the best spectral ratio model is more robust 
than the best single band model with regard to changes in sun angle (season), 
atmospheric transmission, and instrument settings between these two 
overpasses of LANDSAT 7 (July 1 , 2000) and LANDSAT 5 (September 27, 
2000). 

[00143] The same analysis (extracting the best single band 

combination model and the best spectral ratio model) was performed for the 
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LANDS AT 5 data set, which had only 20 samples collected instead of the 30 for 
the LANDSAT 7 overpass. The regression model derived from spectral ratios 
passed the DW test [R 2 (Adjusted) = 63.2%, PC = 1 6.9 + 58.3(R31 ) - 1 08(R42) - 
31 .5(R53) -1 .63(R75)], but was not as strong as the LANDSAT 7 data set model, 
presumably because of the smaller sample size for the L5 overpass data set. 
The single band model was once again undetermined in the DW test. 
[00144] Models for Chlorophyll a derived from both LANDSAT 7 and 

LANDSAT 5 datasets had high R-sqr values, but both of them failed to apply well 
to each other. While there was no reasonable model that could be obtained for 
the IVPF/IVCF ratio from the LANDSAT 7 dataset [R 2 (Adj.) = 35.5%], a strong 
correlation was found between this ratio and a spectral ratio derived from 
LANDSAT 5 dataset [R 2 (Adj.) = 81 .5%]. Application of this model to the 
LANDSAT 7 dataset generated a worse result. The possible reason for the poor 
performance of any model for the IVPF/IVCF ratio applied to the LANDSAT 7 
data set is that the water sample obtained on September 27, 2000 (a LANDSAT 
5 overpass day) was highly dominated by cyanobacteria, as evidenced by the 
much higher IVPF/IVCF (>3.5) than that obtained from the July 1 , 2000 
(LANDSAT 7) data set (<1 .0). 

[00145] The application of the best LANDSAT 7, July 1 , 2000 model for 

PC (using spectral ratios) to the LANDSAT 7 frame of July 1 , 2000 is shown in 
Figure 6, where redder color in the water corresponds to higher amounts of PC 
(in micrograms/liter). Figure 7 shows an image of the same LANDSAT 7, July 1 , 
2000 spectral ratio PC model applied to the LANDSAT 5 frame of Sept. 27, 2000. 
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[00146] Note that there is much higher PC in the September 27, 2000 

image (Figure 7) than in the July 1, 2000 image (Figure 6), which coincides with 
the expected temporal occurrence of blooms of Microcystis spp. that typically 
peak in Lake Erie in late September, for years in which they have been recorded 
to occur. 

[00147] Although all 15 possible non-reciprocal (RJI was considered 

redundant to RIJ) spectral ratios underwent multiple regression analysis, the best 
spectral ratio model in Equation 5 employed only six ratios: R31, R41, R43, R53, 
R73, and R74. LANDSAT TM band 3 is in four of those spectral ratios. Figure 8 
(Green, 2003) shows reflectance spectra of water samples from two lakes, one of 
which contains predominantly chlorophyll a and the other of which contains both 
phycocyanin and chlorophyll a. Note that phycocyanin displays its most 
characteristic reflectance minimum within the wavelength limits of TM bands 2 
and 3, but the gap between the two bands contains more phycocyanin absorption 
than either band 2 or 3. One of the spectral ratios in the PC algorithm (Equation 
5) is R31, which is understandable because the average of both spectral curves 
in TM band 1 (see Figure 8) are approximately equal, yet the reflectance is 
significantly higher for the lake water sample with little or no phycocyanin, 
compared to the one with elevated phycocyanin, in TM band 3. The other ratios 
all involve TM bands 4, 5, and 7, but neither of the spectra in Figure 8 extends to 
those wavelengths. It seems reasonable, however, to assume that these 
reflective infrared bands are "sounding" near-surface water depths, at least in the 
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top half-meter, as the extinction coefficients of water increase with increasing 
wavelength. 

[00148] Figure 9 is a plot of relative phycocyanin content and turbidity 

from the 30 water samples collected on July 1 , 2000, showing that the two 
parameters are somewhat correlated. Because of this, one may search for the 
best spectral ratio model for turbidity, as reported herein from the July 1 , 2000 
water samples that passed the DW test, so that one may compare the turbidity 
and PC patterns in a body of water, as was done with respect to Lake Erie. In 
this example that search resulted in the following equation, with R 2 (Adj.) =85.2% 
and S=1 .579 NTU (Nephlometric Turbidity Units), which is about 9% of the total 
turbidity range: 

[00149] Turbidity (NTU) = -17.2+27.7(R32) 

Eqn. 6 

[001 50] Note from Equation 6 that this turbidity model includes only the 

R32 spectral ratio, which was not included in the PC algorithm of Equation 5. 
Muddy water has higher reflectance in the visible red (TM band 3) than the 
visible green (TM band 2) wavelength regions, so the fact that R32 is positively 
correlated with turbidity is not surprising. Figure 10 shows actual versus 
predicted (Equation 6) turbidity for the July 1, 2000 LANDSAT 7 frame. 
[00151] This turbidity model may then be applied to a small portion of 

data, such as the July 1 , 2000 frame of LANDSAT ETM+ data, which is shown in 
Figure 1 1 with the PC model image of the same area for comparison. The areas 
of highest PC occur in relatively high turbidity zones, but may not correspond 
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with areas of highest turbidity, such as those which are in the Maumee River, 
upstream (SW) of the mouth. This, plus the fact that none of the ratios are 
repeated between Eqns. 5 and 6, are evidence that one may indeed map PC and 
turbidity separately, though the two parameters tend to be somewhat correlated 
in the original water samples. 

[00152] The detection of algal pigmentation signatures from space by 

satellite remote sensing permits the present invention of large areas of aquatic 
systems synoptically. Guided by advances in pigment analysis, the compositions 
of algal populations in aquatic systems are theoretically detectable by remote 
sensing. To date, the application of remote sensing to aquatic systems has 
focused mainly on the detection of chlorophyll a (Abbot and Chelton, 1991 ; Smith 
and Baker, 1982), the pigment common to all phytoplankton. Several satellite- 
borne sensors have been used for this purpose, including Coastal Zone Color 
Scanner (CZCS), Airborne Ocean Color Imager (AOCI), SeaWiFS and Thematic 
Mapper (TM). The advantage of CZCS, AOCI, and SeaWIFS is that they have 
more spectral bands than TM and can probe the "gap" between the visible green 
and red bands (2 and 3) of TM. The advantages of LANDSAT TM over CZCS 
and SeaWiFS is that TM has a spatial resolution of 30 m, which is ideal for lakes 
and coastal zones that are difficult to study with the 1 -kilometer spatial resolution 
of CZCS and SeaWiFS. In addition, TM has longer wavelength bands (5 and 7) 
extending to 2,350 nm wavelength. 

[00153] Several groups have successfully used TM data for detection 

of phytoplankton (Galat and Verdin 1989, Gitelson et al. 1993, Richardson et al. 
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1991 ). LANDSAT has also been used to detect cyanobacterial blooms, albeit it 
on the basis of chl a distributions (Galat et al. 1990). As a result, that particular 
approach is of limited value, due to its inability to discriminate between 
phytoplankton groups. The present study is significant because it represents the 
first successful effort using LANDSAT TM to detect phycocyanin, a pigment 
specific to cyanobacteria as well as some cryptophytes. Until now, it was 
believed that the TM sensor was not suitable for detecting accessory pigments, 
including phycocyanin (with an absorption feature at 620 nm), because the TM 
bandwidths range from 20 nm to 80 nm. There have been several previous 
studies using hyperspectral sensors, such as AVIRIS and CAMS (both airborne 
sensors), to evaluate cyanobacterial populations (Millie et al. 1992; Richardson, 
1996). The present findings have important implications for the future application 
of LANDSAT TM data on assessment and prediction of water quality of aquatic 
systems, not only by primary production estimates in which chlorophyll a serves 
as an important indicator, but also by the mapping of noxious cyanobacteria 
blooms in Lake Erie and other freshwater lakes and tributaries. 
[001 54] Another important finding of this study is that the model of the 

preferred embodiment derived from dark-object-subtracted spectral ratios is 
much more robust than any model one might derive from a combination of single 
spectral bands. The spectral ratio models obtained from both the LANDSAT 7 
and the LANDSAT 5 data sets passed the DW test and could be applied to a 
data set collected at a different time with reasonable accuracy. In contrast, the 
single band models were undetermined in the DW test (meaning that there were 
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autocorrelation problems) and these models could not be used accurately on a 
frame collected on another date. Therefore, the spectral ratio models were more 
robust than the single band combination models. This is a result that should hold 
regardless of the sensor data employed, not just to LANDSAT TM data, because 
it is based on the empirical removal of atmospheric haze prior to ratioing, which 
makes the spectral ratios directly proportional to the actual reflectance ratios of 
whatever the sensor is observing. If the reflectance spectrum is known for one 
area in the scene that does not change in time, such as the Marblehead quarry 
on the southern shore of Lake Erie in our case, the proportionality constants can 
also be determined and each ratio can be normalized for multiplicative changes 
in the ratios. In this experiment, this was not done, but that ratio normalization 
procedure would have improved the results in Figure 4, making the PC model 
even more robust. The algorithm employed by Subramaniam et al (2001) for 
mapping Trichodesmium spp. in the ocean mixed the use of single bands and a 
ratio of two-band differences, which is not likely to be as robust as the dark- 
object-subtracted ratio algorithm produced by the method employed for 
phycocyanin in this investigation. 

[001 55] The phycocyanin model from July 1 , 2000 was then applied to 

a July 16, 2002 LANDSAT 7 frame of the adjacent frame center (P. 19, R. 31) in 
our first attempt at employing the model to detect cyanobacterial blooms. The 
results of that experiment are shown in Figure 12, with the color red 
corresponding to the highest levels of relative phycocyanin content. Note that 
the greatest concentrations of phycocyanin are indicated in Maumee Bay and 



37 



north of Sandusky Bay, both in the western basin of Lake Erie, as well as along 
the northern shore of Lake St. Clair. 

[00156] Figure 13 shows the same type of image for a LANDSAT 7 

overpass on August 1, 2002. In both images, red corresponds to highest 
phycocyanin content (from 4.98-12.00 pg/L), which indicates that phycocyanin 
content in Lake Erie was greatly reduced on August 1, 2002, compared with July 
16, only sixteen days earlier. A synoptic survey of Maumee Bay conducted in the 
research boat on August 1 , 2002, however, indicated the presence of 
Microcystis, but the earlier July 16, 2002 overpass had no corresponding 
validation of the presence of Microcystis or any other cyanophyte. The August 1 , 
2002 sighting of Microcystis colonies floating at the surface amounted to about 
one or two colonies per square meter, a low count. Figure 14 shows PC for an 
August 8, 2002 LANDSAT 7 frame of Path 20, Row 31 , with an optimal stretch of 
the image (red is 11.57-12.51 pg/L, orange is 10.25-11.57 pg/L, yellow is 8.00- 
10.25 pg/L, green is 5.62-8.00 pg/L and blue is 0-5.62 pg/L. If this image had the 
same stretch limits as Figures12 and 13, all colors in the Figure 14 image would 
have been red except for some of the dark blue pixels, indicating far more PC in 
westernmost Lake Erie than on either July 16 or August 1 , 2002. There were no 
LANDSAT frames of P19R31 collected after August 1, 2002 sufficiently cloud- 
free for PC analysis until December, 2002. 

[00157] Local media reported on September 17, 2002 that a substantial 

bloom of Microcystis had enveloped parts of the western basin in the vicinity of 
the Lake Erie islands (Toledo Blade, 2002). Although it is tempting to suggest an 
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emergent Microcystis bloom on July 16, 2002 was documented, the reduction in 
PC (as measured by the July 1 , 2000 PC model) evident on August 1 , 2002 does 
not support this. It is possible that the model was improperly implemented on the 
August 1, 2002 frame, but those results have been checked and re-checked for 
that possibility and the implementation was found to be true. Alternatively, the 
variation in apparent PC levels may simply represent natural fluctuations in 
cyanophyte populations. Recently, Brunberg and Blomqvist (2003) reported 
such fluctuation in a pelagic population of Microcystis in a Swedish lake. The 
variation might also be the result of physical mixing that dispersed the population 
evident on July 16, 2002 throughout the water column to depths inaccessible to 
detection by LANDSAT. Another possible natural explanation is that 
phycocyanin produced by organisms other than cyanobacteria were responsible 
for the higher PC content on July 16, 2002, and the Microcystis bloom may not 
have occurred until August. Four genera of cryptophyte algae produce 
phycocyanin (Hill and Rowan, 1989; Deane etal. 2002), including one, 
Chroomonas norstedtii, that has been reported to occur in the western basin of 
Lake Erie (Makarewicz, 1993). Subsequent microscopic and lab analysis of 
water samples collected near July 16 indicated the presence of cryptophytes, 
although their specific identity was not determined. Thus, it is possible that a 
phycocyanin-containing cryptophyte was at least partially responsible for the 
phycocyanin signature observed on July 16 (Figure 12). Similarly, single-celled 
cyanobacteria (such as Synechoccus sp. and Synechocystis sp.) could have 
provided the basis for the elevated phycocyanin signature, yet would not have 
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been readily observed as part of routine microscopic examination of the water 
samples from that date. 

[001 58] In either case, it seems likely that at least the August 8, 2002 

image (Figure 14) was the beginning of a large Microcystis bloom reported on 
September 17 in local media. However, the occurrence of a reduced PC content 
on August 1 , 2002, presents complications that may not allow one to detect early 
stage cyanobacterial blooms. 

[00159] Since the extensive bloom of Microcystis in Lake Erie in 1995 

(Brittain et al. 2000; Budd et al. 2002), researchers have continued to examine 
various potential bloom-causing conditions such as changes in nutrients, 
temperature, light levels, or even selective feeding - and subsequent rejection as 
pseudofaces by zebra mussels (Vanderploeg et al. 2001). 

[00160] The present invention thus is a method that employs LANDSAT 

TM remote sensing technology for quantitatively mapping relative phycocyanin 
content of freshwater lakes from space or other remote locations. From water 
samples collected during two LANDSAT overpasses, a spectral ratio algorithm 
(Equation 5) for PC with a standard error of 0.592 jjg/L (about 14.8% of the PC 
total range on July 1 , 2000) was created that predicted the PC values on a 
LANDSAT 5 overpass with an error of 2.00 pg/L (about 18.2% of the PC total 
range on September 27, 2000). This error is sufficiently small to permit us to 
map large increases in PC that occur seasonally, and hence, to map 
cyanobacterial blooms, with the assumption that they also create differences in 
PC of similar or larger magnitude. Because the PC and turbidity data from the 
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July 1 , 2000 water samples were somewhat correlated, we also created a 
turbidity model and showed that the highest PC values and the highest turbidity 
values did not coincide, though the highest PC values were found in relatively 
high turbidity areas. This, plus the non-overlap in spectral ratios employed by the 
PC and turbidity algorithms, provided evidence that we were mapping PC and 
turbidity separately. Thus, the concentration and spatial distribution of 
cyanobacteria can indeed be assessed by use of the LANDSAT TM sensor's six 
visible/reflective infrared spectral bands 

[00161] The second question was not completely answered by this 

investigation, though the August 8, 2002 LANDSAT 7 frame for P20R31 
(including the city of Toledo) showed a marked increase in PC (using the best 
July 1 , 2000, PC spectral ratio model) from both July 16 and August 1 , 2002, and 
local media on September 17, 2002 announced a major Microcystis bloom. The 
drop in PC content from July 16 to August 1 , 2002 presents a complication that 
has yet to be resolved, though it could have been caused by phycocyanin from 
other phytoplankton (cryptophytes or other cyanobacteria) on July 16, 2002, or 
as a result of natural fluctuations in a population of Microcyctis. 
[00162] The present invention takes advantage of two technology- 

based approaches: the use of remote sensing, which quantitatively measures 
light reflected from the surface of the earth, as a tool to study regional-scale 
aquatic ecosystem dynamics, and the refinement of techniques to identify and 
quantify algal pigments (Richardson, 1996). 
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[001 63] Accordingly, the findings from practice of the method of the 

present invention are as follows: 

[00164] The LANDSAT TM sensor can be used to evaluate water 

quality by means of detecting accessory pigments of algae, such as 
phycocyanin, despite the relatively wide range of the sensor's spectral band 
width. Detection of phycocyanin may serve as a tool to detect and map 
cyanobacterial blooms, including blooms of the noxious Microcystis. 
[00165] LANDSAT 7 and LANDSAT 5 datasets can be used together to 

enable a more timely assessment and early stage detection of noxious 
cyanobacterial blooms in Lake Erie by providing an 8-day monitoring interval 
instead of the 16-day interval that only one satellite provides. 
[00166] Spectral ratio multiple regression models are more robust and 

reliable than single band multiple regression models for mapping relative 
phycocyanin content from LANDSAT TM data. 

[00167] The 30-meter spatial resolution of LANDSAT TM bands 1-5 

and 7, which were the only ones used in this experiment, is adequate to reach 
upstream into tributaries of Lake Erie for the measurement of PC and turbidity. 
[00168] Mass blooms of cyanobacteria occur in freshwater and 

estuarine ecosystems throughout the world. These results, based on 
cyanobacteria in Lake Erie, show that one may employ the remote sensing 
technology to detect and map potential toxic cyanobacterial blooms elsewhere. It 
is believed that the method of the present invention will allow for the extraction of 
valuable information for assessment, monitoring, and early stage detection of 
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cyanobacterial blooms, water quality, and aquatic ecosystem health from 
LANDSAT TM data on regional and global scales. 

[00169] The present invention thus helps expand the understanding of 

the cyanobacterial bloom-forming process and its temporal/spatial distributions in 
Lake Erie and other systems. Further, one may gain a comprehensive 
understanding of the quantitative connections between algal accessory pigments, 
spectral signatures, and assessment of aquatic ecosystem function by the 
detection of pigments. The present invention may also provide an appropriate 
platform for further development, testing, and improvement of algorithms applied 
to monitor and predict water quality, which are analytically derived and have true 
multi-temporal robustness. 

[00170] The color key for the total phycocyanin assay data shown in 

Figures 2 -9, 14: 



Color 


Colonies (Phvcocvanin) 
Der 100 Ml of Water 


Red 


14000-18300 


Orange 


10800-13900 


Yellow 


10200-10700 


Yellow-Green 


9600-10100 


Green 


7000- 9000 


All Blues 


0 - 6900 
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